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^ ; Abstract 

On . The vacuum expectation values of the so-called Q-operators of certain integrable 

quantum field theories have recently been identified with spectral determinants of 
particular Schrodinger operators. In this paper we extend the correspondence to 
the T-operators, finding that their vacuum expectation values also have an inter- 
Mh. pretation as spectral determinants. As byproducts we give a simple proof of an 

earlier conjecture of ours, proved by another route by Suzuki, and generalise a 
problem in VT symmetric quantum mechanics studied by Bender and Boettcher. 
We also stress that the mapping between Q-operators and Schrodinger equations 

K^ . means that certain problems in integrable quantum field theory are related to the 

^ ' study of Regge poles in non-relativistic potential scattering. 



e-mail: p.e.dorey@durhain.ac.uk 
e-mail: tateo@wins.uva.nl 



1 Introduction and review 

An unexpectedly precise connection between certain two-dimensional quantum field 
theories and the quantum mechanics of anharmonic oscillators was uncovered in ffl . The 
quantum-mechanical side of this correspondence concerned the Schrodinger problem 

H^l^(x) = f-±- + \x\^A^(x)=E^ix) (1.1) 

on the real line x £ (—00,00) , and the associated spectral determinant Dm{)^)^ defined 
for M > 1 as 

DM{E) = DMm\{(l-§-^, (1.2) 

with the product running over the set {E^} of eigenvalues of H. The normalisation 
-C^Af (0) = sin(7r/(2M-|-2))~^ is automatic if Dm is defined as a zeta-regularized func- 
tional determinant (see ||2|), and is also the natural choice for current purposes. However 
we have chosen to negate E as compared to the conventions of [|l|,y] , so that the zeroes of 
Dm{E) as defined here coincide with the eigenvalues of H rather than their negatives. 

The first result of [l| was a mapping between the functional relation satisfied by 
D2{E), discovered in Q, and another set of functional relations, known as a Y-system, 
which had previously arisen in the study of perturbed conformal field theories Q . (The 
precise relation given in [^ had in fact also cropped up in the context of integrable lattice 
models, in 0.) One consequence was a novel expression for D2{E) as the solution of a 
certain non-linear integral equation, a result that was fixed uniquely by the functional 
relations and known analyticity properties on the two sides of the correspondence. The 
move to higher integer values of M brings a considerable increase in the complexity of 
the functional relations. In [|^ the spectral problems for M > 2 were related to certain 
further sets of functional relations, associated with the Y-systems, called T-systems, 
but this result rested on a combination of analytical and numerical evidence without 
being proved. This gap was filled by Suzuki in Q; an alternative proof will be given 
below. 

The other strand in [|^ was a relation between the even and odd spectral subde- 
terminants for the problem ( |1.1| ) and the vacuum expectation values of objects known 
as Q-operators [^,H]! taken at certain values of the 'Virasoro parameter' p. The sub- 
determinants are defined via a split of the eigenvalues according to the parity of their 
eigenf unctions, as 

^m(^)=^m(o) n (i-#) (1-3) 

I even \ >^ / 



' odd 



(so that Dm{E) = D^j{E)D^.^{E)), while Q-operators Q+ and Q_ were introduced, 
in the form that we shall need, in ^. They are continuum analogues of the Q-matrices 
of Baxter |M|. To state the result of [l[ more precisely, we recall that in the construction 
of Q, the operators Q-t act on Virasoro modules labeled by p, with central charge 



c = 1 — 6(/3— /3~^)^, (3 being a further parameter. In each module the 'vacuum' (highest 
weight) state \p) has conformal dimension {p/ 13)'^ + (c— 1)/24, and is an eigenstate of 
the Q-operators. Setting 

^±(A,p) = A^2p/;3^(p|Q±(A)|p), (1.4) 

the eigenvalues A±{\,p) are entire functions of the square of the spectral parameter A, 
and we have 

A±{X,P^/A) = a^Dl,[\^/u^) (1.5) 

where M = /?~^ — 1, v = (^/3^)^~^ r(l— /3^)~^, and the proportionality constants a 
are rederived below. 

In contrast to the first result, M does not have to be an integer for this to hold, and 



so the modulus signs in (IJ.) are in general obligatory when the problem is posed on the 
full real line. However (as exploited in |^] for M half-integral) we can alternatively set 
up the problem on the half-line (0, c«) with the potential x , so long as the boundary 
condition at x = is chosen correctly. The even wavefunctions are picked out by the 
Neumann condition V'2fc(^) ~ '-'' ^'^'^ ^^^ ^^^ ^^ ^^^ Dirichlet condition '(/'2fc+i(0) = 0. 
These two spectral problems thus yield the even and odd spectral subdeterminants for 
the whole-line problem directly. 

The problem of finding spectral determinants related to the Q-operators at general 
values of p was addressed in [10|, where it was found that the problem (|1.1|) should be 



modified to 

-^ + x^^^ + ^)^(x)=i^^(x) (1.6) 

on the half-line (0, co). Imposing the two possible power-like behaviours at a; = 0, 
namely il^{x) oc x^~^^ and ip{x) oc x~\ results in two different spectral problems, and 
their spectral determinants are proportional to the functions j4+(A, (2/-|-l)/9^/4) and 
A_(A, (2/-hl)/?V4) respectively. 

Much of this work rested on the so-called quantum Wronskian relation satisfied by 
the operators Q+ and Q_ |g|. But a key feature of [0,^ was another functional equation, 
relating the Q-operators to further operators T(A), sometimes called quantum transfer 
matrices. Called the T-Q relation, it is 

T(A)Q±(A) = Q±(g-iA) + Q±((?A) , (1.7) 

where 

g = e^"^'. (1.8) 

The vacua \p) are also eigenstates of the T(A), and if we set T{X,p) = {p\T{X)\p) then 
the T-Q relation for these vacuum eigenvalues can be written as 

T{X,p)Ai{X,p) = e^2-M^(g-iA,p) + e±2-Pyl±(gA,p) , (1.9) 

an equation that was also obtained in pl|. 



In this paper we point out that this relation also has a natural interpretation in the 
context of the Schrodinger equation, thus finding a role for the T operators at general 
f3 and p in the 'Schrodinger picture'. This leads to an alternative derivation of the 
results of [0,|l^, and a novel interpretation of the fusion relations and their truncations. 
It also has the bonus, for integer values of M, of providing a simple proof of the T- 
system conjecture alluded to above. All of this material is contained in §§|2|-y, while 
^ applies the results to a problem in quantum mechanics, and §^ discusses duality 
properties. Finally, ® contains our conclusions, and an appendix details the calculation 
of a certain asymptotic used in §^. 

2 The general T-Q relation 



The results of |12| provide a convenient framework for our discussion, and we begin with 
a quick summary of some of this material. 
Consider the differential equation 

-^+P{x)^i;{x)=0 (2.1) 

for general complex values of x. For our purposes we must set P{x) = x — E + 

l{l+l)x~'^ . In [|l^] P[x) was taken to be a polynomial in x, restricting 2M to be a 

positive integer, and /(/+1) to be zero. Allowing 2M to be a general real number larger 

than —2, or /(/+1) to be nonzero, should not change the result quoted below in any 

essential way, though it will usually introduce a branch point in ip{x) at x = 0. Modulo 

this proviso, we have []T^] : 

• The equation ( p.lD has a solution y = y{x, E, I) such that 

(i) y is an entire function of {x, E) (though, for the reason just mentioned, x must in 

general be considered to live on a suitable cover of the punctured complex plane). 

(ii) y and y' = dy/dx admit the asymptotic representations 

y ~ x-*^/2 exp(-^x*^+i) (2.2) 

y' ~ -x*^/2 exp(-^x^'^+i) (2.3) 

as X tends to infinity in any closed sector contained in the sector 

. . Stt , ^, 

argx < . 2.4) 

' ^ ' 2M + 2 ^ ^ 

Furthermore, the solution y is uniquely characterised by this information. (Note though 
that the asymptotic ( |2.2| ) is only valid as given if M > 1. Extra terms appear for 
M < 1, consistent with the WKB result that the solution which decays as j; — > +oo is 
asymptotically proportional to P(x)~^'*^ exp(— J^ P{t)^''^dt) .) 
Let Sk denote the sector 

2kiT vr 



From (2.2) it follows that y tends to zero as x tends to infinity in Sq, and to infinity 
as X tends to infinity in 5_i and Si. More technically, one says that y is sub dominant 
in 5o, and dominant in S±i. To find solutions sub dominant in other sectors, consider 
y{x) = y{ax, E, I) for any constant o. This function solves the equation 

--^ + a-^-^V--a^^ + i(i±H)^(.) = 0. (2.6) 

Thus if a^*^"*"^ = 1, y{ax,a^'^E,l) is another solution to the original problem (|2.1| ). 
Setting uj = exp(7ri/(M+l)) we therefore have a set of solutions 

yk = yk{x, E, I) = uj^'^y{u:-^x, co^^E, I) , (2.7) 

with yk subdominant in Sk and dominant in Sk±i- (Our convention differs from [Q by 
the factor of a;'^'^, which is included for later convenience.) 

A consequence of these facts is that each pair {yk,yk+i} provides a set of linearly 
independent solutions to the second-order equation (p.lj), and any other solution can be 
expanded in terms of them. In particular, 

yk-i{x, E, I) = CkiE, l)yk{x, E, I) + Ck{E, l)yk+i{x, E, I) . (2.8) 

The functions Ck and Ck are called the Stokes multipliers for yk~i with respect to yk and 
yk+i. It follows from 1^) that Ck{E,l) = Ck-i{uJ^E,l) and Ck{E,l) = Ck-i{uJ^E,l). 
For brevity we will write Cq and Cq as C and C respectively. (Again, we differ slightly 



from the conventions of |12|, where the abbreviations C and C were instead reserved 
for Ci and Ci.) 

The Stokes multipliers can be expressed in terms of Wronskians. Recall that the 
Wronskian W[/, 5] of two functions f{x) and g{x) is defined as 

W[f,g] = fg'-fg. (2.9) 



If / and g both solve a Schrodinger equation such as (|2.1|) , then VF[/, 5^] is independent 
of x; furthermore, it vanishes if and only if / and g are linearly dependent. Taking the 
Wronskian of (^^) at A; = with yi and yo shows that 

where we used the abbreviation Wk-^^^k2 ^^^ ^^[yfci) 2/^2]- These Wronskians are entire 
functions of E and I. Since yo and yi are independent, VFo,i never vanishes, and C and 
C are also entire. 

In fact, all of the Ck are identically equal to —1 |l2|. This follows from (|2.1C| ) and 
the relations Wk-^+i,k2+iiEJ) = '^ki^i'^'^^^^) ^^^ Wq^i{E,1) = 2i. (The second of 
these is found by evaluating Wo,i as x tends to infinity in the sectors Sq or Si, where 



the asymptotic behaviours of yo and yi and their derivatives are determined by (2.2) 



and Q.) Since y_i{x,E,l) = yi{x*,E* ,1*)*, it also follows from ( p^ that C{E,l) 
is real whenever E and / are real. 

The basic Stokes relation (12. 81) at /c = is therefore 



C(S, /)yo(x, E, I) = y_i(x, ^, /) + yi(x, ^, (2.11) 

with 

CiE,l) = ^W.i,i{EJ). (2.12) 

2i 



If ( p. 11 ) is rewritten in terms of y it becomes 

C{E, l)y{x, E, I) = uj-^/'^y{uJx, uj-^E, I) + u^/^y{uj-^x, lo^E, I) (2.13) 



With X formally set to zero, this has exactly the form of ( |l.9D for A^{X,p), albeit at the 
specific value /?^/4 of p for which e^'^^^ = lo^I'^ . (It also matches the T-Q relation for A- 
at p = — /3^/4, but since j4_(A,p) = A+(A,— p) P this is not an independent result.) 
However this tactic only works when /(/+1) = 0. Otherwise, the resulting equation in 
E is either trivial or meaningless: if /(/+1) < 0, that is —1 < / < 0, then y{x=0, E, I) is 
identically zero, while if /(/+1) > then y{x=0,E,l) is almost everywhere infinite. 

The problem arises because any solution to ( |2.1| ) is a linear combination of one 
solution, tp^ , behaving near x = as x^"*"^, and one, tp" , behaving there as x~'. Both 
of these are zero at x = if /(/+1) < 0, whilst one or other is infinite if /(/+1) > 0. 
However, rather than considering the functions y{x, E, I) at x=0 directly, we can take a 
hint from the result of |l^ and project onto either ip^ or 7p~ . We choose to fix ip~^ by 
the X ^ asymptotic 

4^+{x,E,l) r.x^+'^ + 0{x^+^). (2.14) 

Since ip~ , the other solution, behaves as x~\ this only determines ip~^ uniquely ii^el > 
—3/2. If necessary, ip~^ can be defined outside this domain by analytic continuation. 



In particular, since / only appears in (2T) in the combination /(/+1) , we can continue 
from / to —l—l and define ip~ as 

'ip-{x,E,l)='ip+{x,E,-l-l). (2.15) 

This procedure does bring some subtleties, to which we shall return §S below, but they do 
not affect the arguments of this section. In discussions of the radial Schrodinger equation 
(see, for example, chapter 2 of [13| or chapter 4 of [|14| ) ■ip^{x,E,l) for "Siel > —1/2 is 



sometimes called the regular solution 



In analogy to (|2.7D, we define 'shifted' solutions tpf. 



i>^ = Vfc (x,^,/) = LJ^/^^j^{uj-''x,uj^''E,l) . (2.16) 

These also solve the original problem (^.l|). By considering the x ^ limit it is easily 
seen that 

iP^{x, E, I) = wT'=('+i/2)V)±(x, E, I) . (2.17) 



We also have W[y^^^^,ilj^^^^]{E,l) = W[y^^,ilj^^]{uj^E,l) , so 

= uj^''^^+^/^^W[y,ij^]{oj^''E,l). (2.18) 

We can now take the Wronskian of both sides of ( ^.lll ) with ip . Defining 

D^{E,l) = W[y{x,E,l),ij^{x,E,l)] (2.19) 

and using equation ( |2.18| ), the Stokes relation ( |2.1lD becomes 

C{E, l)D^{E, I) = w^(i/2+02)T(^-2^^ i^ ^ a;±(i/2+0i^T(^2^^ ^^ ^2.20) 



q2. 



and (|l.9|) has indeed been matched, provided uj is equal q (that is, e^'^/i^^+'^> = e^'^P ^^ 
and uj'^^l'^ is equal to e^'^*^. These are the same relations between M and /3, and / and 
p, as obtained in [Q,|l3li found here by an alternative route. 



To establish the precise relation between the functions appearing in equations (1.9) 
and ( p. 201 ), we can use the fact that, when combined with certain analyticity proper- 
ties, T-Q relations of this kind are extremely restrictive [Q, |,15|. Since D~^{E,l) = 



D~(E, —l—l), we need only consider D~{E, I). In addition to (2.20), we have 

(i) C and D~ are entire functions of E ; 

(ii) If / is real and larger than —1/2, then the zeroes of D^ all lie on the positive real 
axis of the complex- £■ plane; 

(iii) If — 1— M/2 < / < M/2, then the zeroes of C all lie away from the positive real 
axis of the complex- £■ plane; 

(iv) If M > 1 then D~ has the large- -E asymptotic 

logD-{E,l)r^^{-Er , \E\^oo, \aTg{-E)\<7r (2.21) 

where // = (M+1)/2M and 

/•OO -1 

ao = 2 y^ [(t^A^ + 1) V2 _ t^]dt = -^r(-i - ^)r(l + ^) ; (2.22) 

(v) If ^ = then 

1 21+1 I 1 

D-{0,1) = ^r(l + ^^) (2M+2)''^^+'^ . (2.23) 



Property (i) follows from the definition ( 2.19| ) of D as a Wronskian, given that the 



functions involved are themselves entire functions of E. Property (ii) is also straightfor- 
ward, since a zero of D~{E, I) signals the existence of an eigenfunction for ( [L.6| ) at that 
value of E, decaying as x'''^'^ as x ^ 0, and exponentially as x ^ +oo. The self-adjoint 



nature of this problem for / > —1/2 then ensures the reahty of these zeroes. For / > 0, 
the potential is everywhere positive and multiplying ( |1.6D by ijj* (x) and integrating from 
to oo shows that all of the eigenvalues E must also all be positive. For —1/2 < / < 
the centrifugal term in the potential, /(/+l)/x^, is negative but the same style of argu- 



ment can be applied to the transformed equation (A. 3), with the conclusion that the 
eigenvalues are again all positive. Property (iii) is more delicate, and further discussion 
will be postponed until §^, where a partial result will be established. Finally, property 
(iv) follows from a WKB analysis, which is outlined in appendix A, and property (v) 
from a mapping of the problem at £^ = into an exactly-solvable case, given in §^ 
below. 

We now claim that for M > 1 and -1/2 </ < M/2 , the T-Q relation (|2.2q ) and 
properties (i)-(v) characterise the functions C{E,l) and D~{E,l) uniquely. Further- 
more, with the identifications 



P' 



M+1 



P 



21+1 

4M+4 



(2.24) 



the same T-Q relation and the same properties (i)"(v) hold for the functions T{\,p) 
and A-\-{\,p) of ||8[, save for A^ replacing E in (i), the asymptotic in (iv) becoming 

logA+(A,p)~(M+l)r(2^)2%(-Ay, lA^l^oo, |arg(-A2)|<^, (2.25) 

and y4-|_(0,p) being equal to one rather than the value given in (v). The asymptotics 
can be made to agree by setting 



A = vE^'^ , 



(2M+2)-^/2A.r( 1 )-i 



2/x> 



(2.26) 



while (v) is fixed by multiplying j4+(A,p) by D (0,/). By the claimed uniqueness, we 
therefore have 

(2.27) 



^±(A,p)., = «^^^((^)',i-^ 



where a^ = D^{0,2p/P^-l/2)-^, and 



T{X,p) 



cm 



2p 



M=f3- 



M=I3- 

|^^.a((^)M?-i 



W 



M=I3- 



(2.28) 



Continuation from / to —1—1 and the identity A_(A,p) = vl-|_(A, —p) were used to deduce 
the relation between yl_ and D^ in (2.27), while the result ( |2.12|) was used to express 
C{E, I) in terms of W-i^i{E, I) in ( 2.28 ). As will seen in the next section, the expression 



in terms of the Wronskian appears to be the most general. 

It remains to discuss the uniqueness property. Observe first that the asymptotic 
( |2.2lD applies as \E\ — > oo along any ray with avgE ^ 0, since the WKB problem only 
has turning-points if E is on the positive real axis. Since the growth of log D^{E, I) 



in this remaining direction is no greater we conclude that the order of D (E, I) as a 
function of E is equal to (M+1)/2M. (The order of an entire function f{z), if finite, is 
the lower bound of all positive numbers A such that |/(^;)| = 0{e^ ) as |2:| = r — > oo. 
See, for example, chapter 8 of ||l^.) For M > 1, (M+1)/2M is less than 1, and 
Hadamard's factorization theorem implies that D~ can be written as an infinite product 
over its set of zeroes {En} of the following form: 



D-{E,l)=D-{0,l)f[(l-^) . 

n=0 ^ "^ 



(2.29) 



As an aside, we note that this is the final ingredient needed to prove that D~ , and 
hence A^, is proportional to a spectral determinant: from the discussion of property 
(ii) above, the zeroes of D^ coincide with those of the relevant spectral determinant, 
and the residual ambiguity of an entire function with no zeroes is reduced to an overall 
constant by ( 2.29| ). (In fact, the normalisation of D^ given by ( 2.231) coincides with an 
earlier 'natural' normalisation for the spectral determinant when M is an integer and / 
is equal to or 1 [y, so it is reasonable to say that D^ is actually equal to the spectral 
determinant.) 

If M < 1 the order is larger, permitting a more complicated prefactor to the product 
in (|2.29| ), and furthermore the asymptotic density of zeroes necessitates modifications 
to the terms in the product to ensure their convergence (again, see [^). For the most 
part we have been excluding such cases from our discussion, but we remark that there is 
evidence that many of our results continue to hold even after the region M > 1 (called 
the 'semiclassical domain' in ||8[) is left. Some of this evidence was given in [y; some 
more can be found in the next section, and in §^. 

Reverting to M > 1, the product representation ( 2.29| ) allows the steps described 
in [|, 17| to be repeated to show that D~ is determined by a nonlinear integral equation 
of the type introduced in [^, 17|. Introduce the function 



d{E,l) 



CO 



21+1 



D-{uj^E,l) 
D-{uj-^E,l) 



(2.30) 



From ( 2.201 ) and property (i), the points at which d = —1 are exactly the zeroes of C 
and D^ . By properties (ii) and (iii), those on the positive real axis are the zeroes of 
D~ , those on the negative real axis the zeroes of C. (The need for property (iii) was 
not mentioned in [^, but if it fails equation ( 2.32| ) below has to be modified, becoming 
the massless version of one of the equations obtained in ||l^.) The large- i? behaviour 
of d{E, I) follows from {^^f\. 



logdiEJ) 



~ < 



-^ibo{l-e-^''/^'^){EY 
-\ibo{l-e^^l^){EY 



^ < 3.Tg{E) < 2^ - 
-^ < arg(i?) < jj^ 

27r 



27r 
Af+1 



M- 
-27r + 



(2.31) 



M- 



- < arg{E) < 



2-K 

AZ+l 



'this asymptotic is only correct for M > 1. The situation for M < 1 is discussed in appendix A 



where bo = 2cos(2f7)ao- (By {EY we imply e'^''''^^^^\E\^'. Thus the first and third 
asymptotics coincide, as indeed they must since d is a single-valued function of E. In 
the language of |17|, they correspond to the 'second determination'.) 

Now trade E for the new variable 6 = iJ,log{i'^E). Then the arguments of ||8[ can 
be followed to show that the function f{9,l) = logd{i'~'^e^''^,l) solves 

f{e, I) = i7r(/+i) - ii bou-^^'e^ + / ip{9-9') log(l+e^('''')) dO' 

Jcx 

if{9-e') log(l+e-^(^'')) dO' (2.32) 



where the contours Ci and C2 run from —00 to +00, just below and just above the real 
0-axis, and 



jke 



sinhf(^-l)A; dk 



^(^)=/ . 2coshf.:inhKfe 2^ ' ^ = ^- ^'-''^ 



As mentioned above, such nonlinear integral equations first arose in [[151 , p^ . We must 
now assume that the solution to this equation is unique. In this ('massless') context such 
an assumption seems reasonable, though we did not attempt a complete proof. (Note 
though that for the closely-related TBA integral equations, such a uniqueness result can 
fail at nonzero complex values of the mass scale, as can be seen from figure 1 of ||l9|| .) 
A knowledge of / fixes D~ up to an overall constant, given the general information on 
the locations of the zeroes of D~ and C contained in properties (i) and (ii) above. In 
turn, this constant is determined by property (v), and so D^ and hence C have been 
determined. While a more direct approach to the main uniqueness claim would perhaps 
be more satisfactory, the fact, already stressed in [Q], that spectral problems can be 
solved with the aid of a nonlinear integral equation is of independent interest, and will 
be put to use in §0 below. 

3 The harmonic oscillator and the free-fermion point 

The case M = 1 falls just outside the 'semiclassical domain' M > 1, /3^ < 1/2 that 
we have been discussing. However, it corresponds to the three-dimensional harmonic 
oscillator, for which the exact wave-functions are known. As explained in, for example, 
[ pO| ] , the general solution to the radial Schrodinger equation for the harmonic oscillator at 
angular momentum I can be written in terms of the confluent hypergeometric functions 
M{a,b,z) and U{a,b,z) (the notation is as in [^). The correctly- normalised solution, 
subdominant in the sector Sq, turns out to be 

y{x,E,l)=x'+'e---'/^U{^{l+l)-§,l+lx^) (3.1) 

and has the asymptotic behaviour 

y{x, E, I) ~ x-^/^+^/^[l + 0(x-2)] exp(-ix2) . (3.2) 



(Notice the extra factor of x^'^ here compared to the ear her formula ( |2.2| ), reflecting 
the fact that the semiclassical domain has been left.) Taking from [^1| the analytic 
continuation formula 

U{a, b, ze^™) = (i_e-2^^fe") |^^^~^^ M(a, b, z) + e-^'^'^Uia, b, z) (3.3) 



and the Wronskian 



W\U{a, b, z),M{a, b, z)] = ^^"'e' , (3-4) 



the function ^W-i^i{E,l) can be calculated explicitly. The result is 

27r 

*»^"<^-"„.: = r{J+i+f)r(H+f)' (3.5) 



At M=l we have u = l/^/ivf, E = inX"^, and / = 4p-l/2, and so from { ^1^ we find 



^^^'^^ /32=i/2 r (i+2p+7rA2) r (i-2p+7rA2) ' ^^'^^ 

This agrees with the result obtained by the authors of [^ for the vacuum expectation 
value of the operator T at /3^ = 1/2, the so-called free-fermion point, for the particular 
value of their non-universal renormalisation constant C. 



One subtlety has been glossed over here. The modified asymptotic (|3.2|) changes the 
value of Wq^i{E,1) from 2i to 2ie'^^^'^, so that C{E,l) is no longer equal to —1, but 
rather — e~'^ ". Hence ( ^.11 ) becomes 



CiE, l)yo{x, E, I) = y^i(x, E, I) + e-^'^/^yi{x, E, I) (3.7) 

and the would-be T-Q relation (|2.20|) must be replaced by 



e^'^'/'CiE, l)D^{E, I) = ^T-('+i/2)+^/2i)T(^-2^^ ;) + u;^^'('+'/'y^/^D^{u;^E, I) 

(3.8) 
with CO = e'^*'^. This precisely matches the 'renormalised' T-Q relation obtained in [^ 
for the free-fermion point, so long as T{X,p) is identified with e'^^^'^C{E,l) at this 
point rather than C{E,l). But when this T is reexpressed in terms of Wronskians, the 
new prefactor is exactly cancelled by the modification to Wo,i{E,l), with the result 
that the general relation ( 2.2^ ) between T and M^-i^i survives unscathed. Considering 



the asymptotics discussed in appendix A, we expect similar modifications to the T-Q 
relation to occur at M = l/(2m— 1) , /?^ = 1 — l/2m , for any positive integer m. 

Now we can return to property (iii) from the last section, which states that for 
I £ (— 1— M/2,M/2), there are no zeroes of C on the positive-real E axis. We will only 



discuss M > 1. Then a WKB analysis as in [22| can be used to show that at large 
\E\ all zeroes must lie near the negative-real axis, and hence satisfy property (iii). (In 
fact, it seems clear that for M > 1 they lie on the negative axis, but this is harder to 
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prove.) There remains the possibility that some zeroes at smaller values of \E\ lie on 
the positive-real axis, and to exclude this at least for M near 1 we use the fact that the 
zero positions at M=l are determined by (|3.5D, and are 



E = -4n-2±(2/ + l) , n = 0,l,... (3.9) 

Thus for M=l and —3/2 < / < 1/2 the zeroes are bounded away from the positive-real 
axis. By uniform continuity in compact domains, this will also to be true for the low- 
lying zeroes if M is sufficiently close to 1. (Note that continuity alone does not exclude 
the possibility of a zero at large \E\ violating (iii) for M arbitrarily near to 1, which 
is why we also had to invoke the WKB argument.) In fact, it appears that all of the 
zeroes of C actually lie on the negative-real axis for M > 1 and / G (— 1— M/2, M/2). 
At / = 0, this is the main conjecture of |22|; some numerical evidence for the claim at 
other values of / is given in §0. However, the given range for / is certainly maximal: 

from the formula 

2/+1 
C(0,0 = 2cos^^7r, (3.10) 

a consequence of the T-Q relation (|2.20D , it follows that C{E, I) has a zero at £^=0 when 
/ = -l-M/2 and I = M/2. 

Finally in this section we calculate D~{0, OIm • ^^ M=\ , ( p.l| ) can be used with 
the relation C/(a, 2a,x^) = -j=x^~'^°'e^ /'^K^_ii2{x'^ /2) , where Ka is a Bessel function of 
second kind, to see that 

y(x,0,OlM=i= ^^'^'^V2+i/4(2:V2). (3.11) 

For zy > and z ^ 0, K^{z) ~ \T{v) {zjiy |l| , so for ^ > -1/2 we have 

y(x,0,/)U,^i~^2'-i/2r(^+i)x-' asx^O, (3.12) 



TT 



and hence, using (Q, ^-(O, OIm=i = ^2^+3/2 r(l+|+i) . 

The general-M case can be recovered by using the fact that at i?=0 a variable change 
relates a solution at arbitrary M and I to a solution at M=l, but with / replaced by 
/' = l'{l,M). Considering that the normalisation at large x must be in agreement with 
( |2.2[) , the relevant transformation is 

y(x,O,0lM=(M^)^/^X-(^^-^l)/V(M^)^/V^^^^)/^O,iiI-i)U,=, (3.13) 

and repeating the steps already described for M=l, the result quoted in the last section 
is recovered: 

21 + 1 . 1 

D-(0,OIm= («")"' = ^^(l + ^)(2M+2)-^+^ + ^ (3.14) 

To find D~^{0, 1), we continue from / to —l—l : 
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21 + 1 I 1 

Z)+(0,OIm= («+)-' = ^r(l - ^) (2M+2)--^+^"'^ . (3.15) 

Notice that 

D-{0, 1)\mD+{0, l)\,, = (a-a+y = {2l+l)/sm ^vr . (3.16) 

When ^=0, D~ and D^ correspond to half-hne problems with Dirichlet and Neumann 
boundary conditions at the origin respectively, and therefore yield the odd and even 



spectral determinants for the whole- line problem (1J-). (This is why the it signs were 
swapped in our definitions between D^ and A^.) Their product is the full spectral 



determinant, and the result ( 3.16| ) at 1=0 matches the formula for Z)j\,/(0) given in |2^,2 



4 Fusion relations and monodromy 

The T-Q relation is not the only functional equation which arises in the context of the 
T- and Q-operators (see, for example, p, 11, 2^-|2^). Further relations are conveniently 
expressed using the 'fused' T-operators Tj, which can be built up by a process known 
as fusion from the basic operator T. Introducing a half-integer valued index j = 0, 1/2, 
1, . . ., the first set of fusion relations, sometimes called a T-system, reads as follows: 

T,{q-y^X)T,{q'/'X) = 1 + T^+y,{X)T^_y,{X) (4.1) 

where To(A) = 1 and Ti/2('^) = T(A). The fused T's can also be obtained from 

T(A)T,(g^+i/2A) = T,_i/2(9^'+'A) + T,+i/2(g^A) (4.2) 

or 

T(A)T,(g-^-V2A) = T^_^/^iq-^-'X) + T^^y^{q-^X) . (4.3) 

The vacuum states \p) are also eigenstates of these fused T-operators. In this section we 
shall show that the vacuum expectation values Tj(X) = {p\Tj(X)\p) arise naturally in 
the context of the Schrodinger equation ( p.lj ) , leading to a reinterpretation of the fusion 
relations and their truncations in terms of the behaviour of solutions to this equation 
under analytic continuation. 

As remarked earlier, each pair of functions {nm, Vm+i} provides a basis for the space 
of solutions to ( |2.1| ). So far, we have only examined the expansion of yk~i in the basis 
{Hk, Vk+i}, but it is natural to ask about other possibilities. To this end, we extend the 
definition ( |2.8| ) of Ck and C^ by setting 

Vk-l = CI 'Vk+m-l + CI 'Vk+m (4.4) 

(so that C^, = Ck and Cj^ = Ck = —1)- The change from the {yk+m-iiUk+m} basis 
to the {yk-i^Vk} basis is then effected by a 2 x 2 matrix C^ as 

yk-l\ _ f^{m) ( yk+m-l\ c'im) _ f C'fc"' C'fe'" A / . .N 

yj-"^' [vk+m) ' "^^^ ~[ci-'^ c^-'^J- ^'-'^ 
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This matrix depends on E and /, but not x. The following properties are immediate: 

CM(E,0 = CK(^=E.0, (4.6) 

Cr^(j ;) . cf'^CJ' -). (.7) 

Further relations reflect the fact that the change from the basis {yk+m+n-i,yk+m+n} 
to {yk+m-i,yk+m}, followed by the change from {yk+m-i,yk+m} to {yk-i,yk}, has the 
same effect as accomplishing the overall change in one go: 

Cf^cg^ = Ci™-^") . (4.8) 

(These express the consistency of the analytic continuations, and can be thought of as 
monodromy relations.) Consider first the case m = 1. We have 






1 ^JW-'^ ct",'') Vca cti. 



(4.9) 



which gives two non-trivial relations: 



■'fc^fc+l ^k+2 ^k 

and 



C.^a - Ct.'^ = Ct'^ (4.10) 



^.^a-^ir2'^=^r'^- (^-n) 



In addition, we have the initial conditions 

■.(0) _ 1 r^(l) 



Cr = l , q'^ = C. (4.12) 

Cf)=0 , C(^^ = -1 (4.13) 



The n = 1 case of (4.11) shows that C\, = — C^ = —C)^ > ^^^ then the more general 



equality 

(5f ) = -Ci"-^) (4.14) 

follows on comparing ( [4.11| ) with ( |4.10[ ) . If we now set 

C^^){E) = C7j")(u;^"+iE) , (4.15) 



then ( [4.10 ) is equivalent to 

C{E)C^'^\uj'^+^E) = C(""1)(c^"+2e) + C("+i)(aj"E) , (4.16) 

and this matches the fusion relation (|4.2| ). Since C^^'{E) = 1 = Tq{E) and, from the 
last section, C^^\E) = C{E) = Ti/2{j^E^''^), this establishes the basic equality 

C^^\e)=T^/,{uE'/'). (4.17) 
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It is easy to check that the fusion relation ([4.3|) emerges in a similar manner from (4.8) 



at n = 1. To recover the T-system (4J), one more piece of information is needed. 
Taking Wronskians in (|4.4|) yields 



^k ~ 7r-^k-l,k+m , C'fc ~ —17:^k-l,k+m-l ■ (4-18) 



An immediate consequence is the recovery of ( [4.14| ), but we also obtain 

-,{-m-2; 
'k+m+1 



ct^ = -<:;? • (4.19) 



Using this result, the n = —m case of ([4.8|), namely C^™' C^,^*^ = 1, implies that 

(j{m-i) (^-i^)C7('n-i) (^^) _ C7(m) (^)c(™-2) (^) = 1 . (4.20) 

Given the identification (4.21), this is the T-system ( |4.1| ), evaluated on the vacuum state 



Finally, the formula ( 4.18 ) allows the function T„/2 to be expressed alternatively in 
terms of a Wronskian: 

Tn/2{^E'^^) = C^'^He) = ^W.,,n{u^-^^'E) . (4.21) 

This will be relevant in §^ below. 

The monodromy relations used so far have all been 'local', in that they can be built 
up from continuations of the functions yk from one sector Sn to its neighbours Sn±i- 
But there is also the possibility to continue all the way round the origin a number of 
times. If the cover of the punctured x-plane on which the yk live closes in a suitable 
sense, this leads to new relations, and these turn out to correspond to the truncations 
of fusion hierarchies in the integrable quantum field theory. 

As an illustration, we will discuss the simplest class of examples, which arise when 
2M is an integer and l{l+l) is zero. Then all solutions to (|2.1| ), and in particular the 
yi^{x,E,l), are single valued on the once-punctured x-plane, and the sector Sn+2M+2 
coincides with the sector 5„. Both y„ and yn+2M+2 ^'^^ subdominant in this sector, and 
so they must be proportional to each other. To find the precise relationship, consider 
their behaviour as \x\ — > cxd in Sn, setting x = pe^'^^'^ ~^^' and letting p -^ oo. From 
(pt) and (13), 

y„ ~ u-'^p'^^^l^ eM-jk-y^^') ' yn^2M^2 - -^"/V-^^/^ eM-jk-y^') (4.22) 
and so 2/„4.2Af+2(^' -^' = ~yni^^E,l). From ( [4.21 ) we immediately deduce 



C(2J^'^)(E) = 1 , C(2A'^+i)(^) = (4.23) 

and hence the relation ( [4.20| ) can be rewritten as 



2M~1 , 



C(™)(u;-^E)C(™)(u;^) = 1+ J] (c(")(^)j (4.24) 



n=l 
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where /„„ is the incidence matrix of the A2M-1 Dynkin diagram. This is the simplest 
example of truncation, and ultimately [^] leads, via the Y functions discussed briefly 
in §^ below, to alternative sets of integral equations of the 'TBA' type [29|. For integer 



values of M, these equations were applied in ||1| to the computation of the energy levels 
of anharmonic oscillators with analytic potentials x . 

At arbitrary rational values of M+1, a similar collapse of the set of solutions 
{yk}kL-oo occurs, but this time on some finite cover of the punctured j;-plane. Again, 
a truncation of the fusion hierarchy results. For nonzero values of /(/+!) the situation 
is a little more subtle, as there is a multivaluedness in the solutions induced by the 
singularity at x = 0. To handle this behaviour the alternative basis {'0~,V'^} is rnore 
appropriate. We leave further discussion of this point to future work, but we expect 



that all of the more general truncations of the fusion hierarchy (4.20) will eventually 
find a geometrical interpretation in terms of the behaviours of solutions to the basic 
differential equation (|2.1|). 



5 The (fused) quantum Wronskians 

The next task is to relate the fused Stokes multipliers C^™"' to the functions D^ discussed 
in ^ Recall that the Wronskian is bilinear in the space of differentiable functions, so 
that given four arbitrary functions f{x), g{x), h{x) and l{x), and arbitrary constants 
a, /3, 7, 5, we have 

W[af + Pg,^h + 6l]=a-fW[f,h] + a6W[f,l]+p-fW[g,h]+p5W[g,l] . (5.1) 

For almost all I (exceptions will be discussed below), the functions {'0~,V'^} intro- 
duced in ( 2.14 ) and ( |2.15| ) provide an alternative basis for the space of solutions to 



the differential equation ( p.l[ ). In particular, using the results VF[y, V'^] = D^ and 

W['il^~,ip+] = 21+1, we have 

{2l+l)y{x, E, I) = D-{E, l)%b'{x, E, I) - D+{E, l)^p+{x, E, I) , (5.2) 

More generally, the shifted solutions defined by ( |2.7| ) and ( |2.16p satisfy 

{2l+l)y^{x,E,l)=D-{uj'^''E,l)ij-{x,E,l)-D-^{u;^^E,l)ij^{x,E,l). (5.3) 

Taking the Wronskian ( ^.3D at k = —1 with the same equation at k = n, shifting E to 
u}^~'"'E and then using the formula ( 4.21| ) for C'"'(£'), property ( ^.1| ) and the results 



W[<P;,^l^t] = (2/+l)u;(^-P)('+V2) (5^4) 

(valid at arbitrary 'shifts' p and /c), we find 

(4/+2)iC(")(E) = Cj("+l)('+l/2)^-(^n+l^^^)^+^^-n-l_g^^) 

-^-i:'^+m+y^)D-{uj-'-'-^E,l)D+{uj''+'^E,l) . (5.5) 
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In the context of integrable quantum field theory, a corresponding set of relations was 
given in j^: 

- q-^^^+^^Plf"A+{q-i-^/^\p)A^{qi+^'^\p) (5.6) 

where j = 0, 1/2, 1 . . . . With the identifications ( p.27| ) and ( |2.28[ ), and using the result 
(|3.16D , the two sets of relations agree if, as before, q = u> and 2p/ fi"^ = I + 1/2. At j = 0, 
( |5.6D is the 'quantum Wronskian' relation of |^, while the n = case of ( |5.5| ) was first 
found (for M an integer and / = 0) in |y]. The match between these two was a key 
ingredient in [|,1C]. 



The 'T-system conjecture' of [|1[, proved in |6[, is a simple corollary of this result. If 
/ = and M is an integer, then (^]^) at n = M becomes 

C7W(i^) = D+{-E)D-{-E) = Dm{-E) , (5.7) 

the second equality following because at / = 0, D^ and D~ are the even and odd 
spectral subdeterminants respectively for the full- line problem ( [1.11 ). Since C^^'^'{E) = 
Tj\,f/2(^-E^^ )) this establishes the conjectured relation between Dm and the vacuum 
expectation value Tjm/2- (Note though the small differences in notation from [|^: as well 
as the negation of E already mentioned in the introduction, a half-integer j has been 
used to index the T operators in this paper, in line with the conventions of |^, 27, 28[ , 
while in |j|] the index was integer-valued. Thus in ||l| the correspondence was with Tm 
rather than T/vf/2.) 

Next we return to the fact that the pair of functions {V'~,V'^} does not always 
furnish a basis for the space of solutions to the differential equation ( p.lj) . The most 
obvious counterexample to such a claim is the point I = —1/2, at which ip~ = ip^. 
More generally, since ip~ is initially defined by analytic continuation in /, there may be 
points at which poles are encountered. If these poles are removed by multiplying by a 
regularising factor, the resulting solution ijj~ may fail to be independent of V'^ at those 
values of I where previously there had been divergences. This is related to the remark 
in Q that at certain values of p the functions A^{X,p) and A-{X,p) may coincide. In 
the context of the radial Schrodinger equation, a discussion of the issue can be found 



in, for example, chapter 4 of [14|. The most direct way to locate the 'problem' values 



of I is probably the iterative construction of ^ and ip~^ given in |3(|. For the potential 



x^^, for which the function U{x) of |30| is equal to x^*^ — E, poles ip~^{x, E, I) occur at 
/ + 1/2 = —mi — (M+l)m2, with mi and m2 non-negative integers. Hence {'ip~,'ip~^} 
fails as a basis at the points 

1 + 1/2 = ±{mi + (M+l)m2) , mi, m2 > . (5.8) 

For integer values of 2M, this is just a standard phenomenon in the Frobenius method, 
which predicts that one of the pair {'ip~ ,'ip~^} may have a logarithmic component when- 
ever the two solutions to the indicial equation a{a—l) = /(/+1) differ by an integer, or 
an even integer when 2M is even. 
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Rather than giving a complete treatment, we wiU focus here on the / = —1/2 case. 
This corresponds to the ground state of the untwisted sine-Gordon model, and so is 
particularly interesting from the field-theoretic point of view. The emergence of loga- 
rithmic corrections can be understood by taking an appropriate limit I — > —1/2. Using 
( |5.2[ ), valid away from the points (^I^), set / = —1/2 -|- e with e tending to zero. Near 
X = we have 



.1^.^^ n-^R _1V . 1/2 ^!±^'_2)-(i^,_ 1)^1/2 



y{x,E,-^+ e) ^ D-iE,-^y x^/^ ^^—;^ - D-{E,-\)x^l^^-^^ (5.9) 



where the primes denote differentiation with respect to / and we used the facts that 
D-{E,l) and D-{E,iy are nonsingular at / = -1/2 to expand D'^iE, -1/2+e) ~ 
D~{E, —1/2) ±D~{E, — l/2)'e. Expanding x"^^ = exp(±elogx) to first order in e, ( |5l9D 
becomes 

y{x,E,-^) r^ D-{E,-^yx^/^ -D-{E,-^) x^/^ log X. (5.10) 

We see that y{x, E, —1/2) has an x^'^ log x component at small x, as expected from the 
Frobenius method. The basis {'ip~,'ilj^} should therefore be replaced at / = —1/2, and 
a suitable choice is the pair {x~ , X^}- 

X+{x,E)r^x^/^ + 0{x^/^) , 

X~ (x, E) ~ x^/2 log X + 0(x^/2) . (5.11) 

The 'Jost functions' ( ^.191) become 

D^{E) = W[y{x, E, -i), x^(x, E)] , (5.12) 

so that D~iE) = D-{E, -i) , and D^{E) = D-{E, -i)' . It is interesting to see the 
effect that these changes have on the basic functional relations. We set Xki^^-^) ~ 
uj^''^X^{'^~^^i^'^^E) s-iid use the Wronskians 

w[x;,xt] = 0, 

W[x~p,X'q] = (9-p)log^, 

W[xt,Xq] = 1, (5.13) 

to find, with D]: = D^{uj'^^E), 

W[yp,yg] = {q-p)D-D~ logo; + {b-D+ - D-D;) (5.14) 

and feeding this into the formula ( 4.21 ) for C'"^ yields 



+ 1 ('^-(cj-"-i^)^+(w"+1e) - D-{io''+'^E)D+{u-''-'^E)') . 
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Note also, from ( |3l^ ) and ( p^ , that 



/ 2 A/f -1-2 2 

^"(0) = \/ ^- , ^+(0) = ^(2+2M)-i/2(log(2+2M) -7s) (5.16) 



where 7^; = — r'(l) = 0.57721... is the Euler-Mascheroni constant. Finally from ( |2.1l| ), 
( |5l0D and ([sTID we obtain 

C{E)b- (E) = L)- (uj^'^E) + L)- {uj'^E) 

C{E)D+{E) = D+{uj-^E) + D+{lo^E) + (d-{uj~^E) - D-{io^E)) logcj . (5.17) 



Thus at / = —1/2 one of the two T-Q relations has to be changed. Similarly, comparing 
( |5.15D with the fused quantum Wronskian ( |5.5| ) shows that the relationship between 
C^'^'{E) and D^{E, I) undergoes a nontrivial modification. The fact that for the ground 
state energy of the untwisted sine-Gordon model things need to be slightly modified was 



already observed, from a completely different angle, in [31|. There, the general solution 
of the sine-Gordon K-system was described in terms of a single quasi-periodic function 
h{E), satisfying h{Eoj'^^'^^'^) = p+qh{E), where p and q can, for the purposes of |31| and 



the example investigated here, be taken to be constants. In order to have the correct 
UV-limit for the untwisted sine-Gordon model (corresponding to / = —1/2), a non- 
vanishing p was required, leading to a log(ii') component in h{E) (see ( |6.14| ) below). In 
a slightly different language, this turns out to be what we have derived here. A match 



with the results of |31| is the topic of the next section, where, as simple byproduct of 



the mapping, a curious property of the spectral determinants is pointed out. 

6 Y-systems and dilogarithm identities 

As mentioned in §Q, there is a second set of functional relations, the so-called Y-system, 
closely related to the T-system discussed in the previous sections. The relation between 
these two systems is 

Yn{E) = C7"+H^)C"-i(£;) (6.1) 

and the y's fulfil the relation 

Yniu;E)Yniu;~^E) = (1 + y„+i(^))(l + yn-i(^)) (6.2) 

For M integer or half-integer and / = 0, this system truncates {Yq{E) = Y2MiE) = 0), 
and it coincides with the ^2M-i-i'elated Y-system discussed in Q. 

On the other hand, the Y-functions are related to the solutions of TBA equations. 
In this framework they encode finite-size effects in integrable quantum field theories, 
and, through the consideration of ultraviolet limits, lead to certain remarkable sum rules 
for the Rogers dilogarithm function involving the stationary {E = 0) solutions of the 
system ( |6.2[ ). For M integer and / = 0, for example, the relevant sum rule is 

fi 2M — 1 / 1 \ 9 A/f 1 

n=l \ ' / 
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where cuv is the central charge of the Z2M parafermionic conformal field theory, L(x) 
is the Rogers dilogarithm 

L(x) = -irdJi^^ + i^^Miz^l , (6.4) 



and the values of the constants Yn{0) involved in ( |6.3D are 

sin(7r ^^) sin(7rT 
^i^ (.2M+2) 



Yn{o)= ' "i:^^; ^ ,' ■ (6.5) 

^ 2M+2 J 



With some additional complications, sum rules similar to (|6.3| ) can be written for any ra- 
tional M [32] and arbitrary I. In ||3^ a generalisation of ( |6.3D involving the -E-dependent 



y-functions was proposed. For an arbitrary solution Yn{E) to (|6.2| ) (but again with M 
an integer and 1 = 0) the result is 

„ 2M-12M+1 . ^ X 

Dilogarithms also appear in certain volume calculations in three-dimensional mani- 
folds (see for example |3j]), and related to this idea is the fact [^ that the general 
solution to ( |6.2| ) can be expressed using cross-ratios 

/ , ,N (a — c)(b — d) , „, 

Surprisingly, the points involved in the cross-ratio can be expressed in terms of a single 
quasi-periodic function of £■ : if h{E) satisfies 

h{Euj^^'^+^) =p + qh{E) (6.8) 

then 

Yn{E) = -{h{u'^E), h{uj-"-^E),h{io-''E),h{uj''+^E)) (6.9) 

solves ( |6.2[) . 

Now we would like to show that the function h{E) is a natural object in the context 
of the Schrodinger equation. Since C^^' (E) = 1 for all E, we can rewrite C^^'^^' (E) as 
C^^^^'[E)/C^^\uj'^~^^E) and, for generic values of I, expand out numerator and denom- 
inator using ( |5.5|) to find 



^ ' D+{io^E,l) {k{uj^+^E)-k{iu^E)) ^ ' 

where 

k{E) = E^+^/^D-{E, l)/D+{E, I) . (6.11) 
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Similarly, 



C(^-^){E) = ^-(m/2)(n+i) D^juJ^'EJ) 



(A:(w"S) - kiuj-'^E)) 



and so, from (|6.1| 



i^n(£^) 



L>+(a;-"-2^,0 (/c(a;-"£;) - k{uj-''^-'^E)) 



k{uj''+'^E) - k{uj-''-'^E) kioj-^'E) - k{uj'^E) 



k{Lo"+'^E) - kiu'^E) k{uj-'^E) - k{ijj-^-'^E) 



(6.12) 



(6.13) 



This was for generic values of I: the story changes at those points where {'0~, V'"'"} foils 
to be a basis. The first example is / = —1/2, where (|5.5| ) is replaced by ( 5.15|). It 
turns out that the final result ( |6.13| ) continues to hold, so long as the definition ( |6.11 ) 
is swapped for 

k{E) =\ogE + 2b+{E)/t)-{E) . (6.14) 



Comparing (6.13) with ( |6.9| ), the function k{E) as defined in (6.11) gives a realisation 
of h{E) for p = and q = e27ri((+i/2)^ while the definition ( |6.14 ) coincides with the case 
q = 1, p = 27ri. Note also that when I = 0, k{E) is essentially the (reciprocal of) the 



alternating, or skew, spectral determinant discussed in 1 23 , 24 1 . It is interesting that the 
result ( |6.6D (and its generalisations) appears to furnish a novel form of sum rule for the 
spectral determinants of Schrodinger equations. 



The T operators, spectral determinants, and a generali- 
sation of a problem considered by Bender and Boettcher 



The result 



^n/2l 



.^1/2) = C7W(^) = 1vF_i„(l<j-"+^^) 
2i 



(7.1) 



from §^ implies that the vacuum expectation values of the T-operators also have an 
interpretation as spectral determinants. The right hand side of (7J) vanishes if and 
only if E is such that the functions y_i and yn are linearly dependent, which in turn is 



true if and only if there is a nontrivial solution to (2J) which decays to zero as x tends to 
infinity in the sectors 5_i and 5„. This is an eigenvalue problem, and the argument just 
given shows that the zeroes of C^""' coincide with the zeroes of its spectral determinant. 
From ( |5.5| ) and (2.21), the order of C^"-* is less than 1 for M > 1, so the discussion around 
( |2.29D for D~ applies equally to C^-^\ and shows that this function is actually equal to 
the spectral determinant, with a normalisation which again turns out to coincide with 
the natural one for 2M an integer and / equal to or 1. (Incidentally, for M integer 
and I zero this gives us another easy proof of the T-system conjecture, obtained more 
indirectly via the subdeterminants D^ in ^ above.) General eigenproblems of this sort, 
where the boundary conditions are specified as \x\ — > oo in two sectors in the complex x 
plane, are discussed in chapter 6 of |]l2| and also in [^, 36|. They are the other natural 
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set of spectral problems associated with (|2.1|), and it is pleasing that they correspond 
so neatly with the fused T-operators. 

As an application, we consider and generalise a problem treated by Bender and 
Boettcher in |^2|. These authors discussed the spectrum of the Hamiltonian p^ — (ix)^ 
with A^ real, corresponding to the eigenvalue equation 

-ip"{x)-{ix)^iP{x) = Ei:{x). (7.2) 

For quantised energy levels the boundary conditions must be chosen appropriately. The 
standard requirement is for ^{x) to tend to zero as |x| — > oo on the real axis, but for 
arbitrary A^ the problem needs to be continued into the complex j;-plane. The real 2;-axis 
is replaced by some contour, and the boundary conditions are imposed at the two ends of 
this contour. The regions where 'ijj{x) can vanish exponentially as |x| — > oo are wedges, 
bounded by what Bender and Boettcher referred to as Stokes lines, with the vanishing 
being most rapid at the centers of the wedges, called by them anti-Stokes lines. (Beware 
that many authors use exactly the opposite terminology - see, for example, [^].) There 



are thus a variety of eigenvalue problems associated with (7^), depending on the choice 
of asymptotic directions for the contour relative to the wedges. The problem studied 
in |2^ was the analytic continuation of the usual harmonic oscillator, so the wedges 
were picked to allow the quantisation contour to run along the real axis when N = 2. 
This selects the wedges centred on the directions 

A^-2 7r , ^ A^-2 7r ,„ , 

c'lpft = —71" + T^ and C^rieht = " TT i 7.3) 

N+2 2 ^ A^+2 2 ^ ^ 

each of opening angle A = 27r/(A^+2). Bender and Boettcher studied the problem both 
numerically and analytically, and found strong evidence for an entirely real spectrum 
for N > 2, with a transition at A^ = 2 below which infinitely-many eigenvalues become 
complex. This they associated with a spontaneous breaking of VT symmetry in the 
quantum-mechanical problem. 

If we set tlj{x) = 'tp{x/i), ( ^ ) becomes 

-i;"{x) + x^ilj{x) = -EiP{x). (7.4) 



and matches our original Schrodinger equation ( p.l| ) with M = N/2, I = 0, and E re- 
placed by —E. Furthermore, it is easily checked that the wedges defined by ( |7.3D become 
the sectors 5-i and Si. Referring back to the discussion following ( [7.1[ ), this means that 
the eigenvalues of the problem posed by (^) and (|7.3|) occur at the negated zeroes of 
C^^'{E), and are thus encoded in the pattern of zeroes of the vacuum expectation value 
of the fundamental quantum transfer matrix T(A) = Tx/2(A). 

In turn, these can be found via the non- linear integral equation ( 2.32| ). The results 



are displayed in figure 1, and agree perfectly with the results of [22|, which were obtained 
by a direct treatment of the differential equation in the complex plane. Note that the 
transition at A^ = 2 corresponds to the point at which the associated sine-Gordon model 
moves from the attractive to the repulsive regime, and that the match with the results 
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Figure 1: Eigenvalues of the Hamiltonian p^ — (ix)^ plotted as a function of N, found via the 
solution of the nonlinear integral equation (2.32|) at I — 0. 



of Bender and Boettcher continues to hold even after this point has been passed, and 
the semiclassical domain N > 2 has been left. 

From the point of view of equation ( ^.32 ), the restriction to / = is unnecessary. 
Once it is relaxed, we are led to the following generalisation of the Bender-Boettcher 
problem: 

- V"(x) - {ix)^^{x) + l{l+l)x~^ = Eij{x) . (7.5) 

This seemingly-innocent modification to the Hamiltonian turns out to have a remark- 
able effect on its spectrum. In figure 2 we show some results from an initial numerical 
investigation. A few features merit immediate comment. At A^ = 2, the exact eigenval- 
ues can always be found by negating formula (^.S|). The behaviour for N < 2 depends 
significantly on the sign of /. The situation at / = —0.025 is shown in figure 2a. It 
might appear similar to that at / = 0, but the connectivity of the states on the plots is 
exactly reversed: the ground state is now joined to the first excited state, the second 
to the third, and so on. As / increases, the level joined to the ground state fiicks up 
through the rest of the spectrum, leaving a restored connectivity pattern in its wake. 
The mechanism should be clear enough from figures 2b to 2d. By the time / reaches 
zero the earlier situation has been recovered, with the ground state no longer connected 
to one of the excited states, but rather diverging to +oo along the line N = 1. Then for 
/ > the ground state eigenvalue ventures for the first time into the domain N £ (0, 1]. 
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2a) I = -0.025 



2b) I = -0.005 





2c) I = -0.0015 



2d) I = -0.001 





2e) I = 0.001 



2f) i = 0.05 



Figure 2: Eigenvalues of the Hamiltonian p^ — (ix)^ + l{l+l)x ^ found via the nonhnear integral 



equation (2.32), plotted as a function of N at various nonzero values of I. 
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(To make sense of the eigenvalue problem in this regime requires its consideration on a 
cover of the once-punctured complex plane, with the sectors 5_i and Si overlapping on 
different sheets.) Figures 2e and 2f show this behaviour. Observe from figure 2f that the 
ground state eigenvalue is zero at / = 0.05, A^ = 0.2. More generally, ( S.lOj) implies that 



the problem has a zero eigenvalue when N = 41. This means that the steeply-climbing 
line on figure 2e should also return to zero, at A^ = 0.004. While our program currently 
breaks down at such small values of A'^, plots made at intermediate values of / (between 
0.001 and 0.05) offer clear support for this scenario. 

For the values of / pictured, the spectrum is entirely real for N >2, thus generalising 
the reality property observed by Bender and Boettcher for / = 0. 

In a subsequent paper Q| further eigenvalue problems were discussed, based on 
Hamiltonians of the form p^ -\- x'^ {ixY , with the two 'quantisation wedges' chosen so as 
to be centred on the positive and negative real axes when e = 0. In these conventions, 
the earlier problem corresponds to -ftT = 1, e = N—2. It is straightforwardly checked 
that, just as T1/2 or C'-^^ encodes the spectrum of the problem for i^ = 1, so the functions 
r^/2 or C^^' encode the spectra of the more general problems. This means that these 
problems, treated separately in [^, are in fact intimately linked together, since their 
spectral determinants participate in the fusion hierarchy, or T-system, discussed in §B 
above. On the other hand, the numerical and analytical work in [^] immediately 
gives us a rather detailed picture of the behaviour of the zeroes of the general fused T 
operators. Worth noting in this regard is the fact that the first transition to complex 
eigenvalues for the Schrodinger problem, that is to complex zero positions for the fused 
T operator, always occurs at e = 0. In the earlier notation of this paper, this point 
corresponds to 2M = N = 2K, and for i^ > 1, it lies well inside the semiclassical 
domain M > 1, /?2 = 1/(M+1) < 1/2. 

As this paper was being finished, a further article by Bender et al appeared ||3^ . 
There a form of 'classical limit' (/3^ — > 0) is discussed. More precisely, the limit is 
M ^ 00 with / = and £ = E/M'^ held fixed. (This differs from the limit treated 
in appendix B of P] where M ^ 00 with I and E is held fixed, and the functions A± 
and Tj, j = 0{M), are treated. Taking the limit with E/M"^ fixed instead of E is 
appropriate for capturing the behaviour of the functions Tj with j remaining finite.) 
For the initial problem, corresponding to K = 1 and the function T1/2 or C'^\ Bender 
et al found 

Ek{M)^{k + \fM'^ , A; = 0,1..., M^oo. (7.6) 



(Note that in [|38[ M is instead used to label the different eigenproblems, the role 
played by K here.) The result ( [7.6| ) suggests that the limiting form for the spectral 
determinant C^'^^S^l) at /=0 is 2cos{'it\J—8) . (The minus sign appears for the same 
reason as before, namely the variable change in going from (|7.2| ) to ( [7. 4] ).) Assuming 



that the fusion hierarchy (4.20) holds in this limit, from it we deduce 



^(i.)(^^0)^^mM+lWE^_ (7.7) 
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This correctly reproduces the spectra for the higher problems that in |38| had to be 
obtained by a series of independent calculations. Finally, we remark that the angular 
momentum term l{l+l)x~'^ drops out of ( |7.5| ) in the limit considered in [38|, and so the 
results (|7.6| ) and ( [7.7D should hold unchanged for / 7^ 0. 

8 Singular potentials and duality 

Up to now the parameter M has been restricted to the range (0, 00), but simple trans- 
formation properties of the Schrodinger equation allow a wider range of M to be accom- 
modated. The idea is to define a mapping between the eigenvalue problem for M > 
and the problem with — 1 < M < 0. Our discussion here is not meant to be particularly 
original, essentially following [^] , but we wish to underline the fact that the relevance 
of integral equations to Schrodinger problems, stressed in M , actually applies to wider 



range of problems than might initially be thought. The first step is to make Langer's |3S] 
variable change 

x = e^ , y{x,E,l)=e'/^'il^{z,E,l). (8.1) 



Then 






+ e2^(*^+^) - Ee'^A ^P{z, E, I) = k^iIj{z, E, I) (8.2) 

with K = —i{l + ^). Now we interchange the role of the two exponential terms, with the 
transformation 

--MTT- + l°gW^ (8-3) 



to obtain 



- -^ - e2^/(^-^+^) + Ee^' ]tP{z,E,l) = k^i;{z,E,l), (8.4) 



where ft = k/{M+1) , I = {l-M/2)/{M+l) and 



E = ^— ^ ^-*'^-i . (8.5) 



It is easy to check that for M = 1 (|8.5D gives an exact mapping between the energy 
levels of the harmonic oscillator at angular momentum / and those of the Coulomb 
potential at angular momentum / = 1/2 — 1/4. Indeed it is straightforward to see that 
eigenfunctions of the original problem are transformed into eigenfunctions for the 'dual' 
problem: reversing the original variable change, we have that 

y{x, E, I) = (Af+l)-i/2^i/^x5P^y((M+l)^'i/2xi/(^^+i), S, /) (8.6) 

solves 

- ^'"^ + ^^ - e\ y{x, E, = 0, (8.7) 
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with M = (M+1) -^ — 1. Furthermore, if y decays as x ^ +cxd then so does y, and if 
y ~ 2;'+^ near x = then y ~ x'"^^. Hence eigenfunctions are mapped to eigenfunctions, 
and the relation ( |8.5| ) is true at the level of the eigenvalues too. A curious consequence 
of this transformation, applied to the problem with M = 1/2 and I = 0, is that the 
solution of the x~^'^ potential with I = —1/6 can be written in terms of the Airy 
function. 

The mapping to singular potentials, decaying (albeit algebraically) at infinity, high- 
lights the fact that the motion of the zeroes of the Q-operators as a function of the 
'twist' p corresponds to the motion of Regge poles for the radial Schrodinger equation 



as the angular momentum / is continuously varied |13,14,4C|. Quite apart from the 
intrinsic curiosity of this fact, it means that the general results of, for example, chapters 

8 and 13 of [0| can be applied to deduce features of the motion of the zeroes of the 
Q-operators as p and /?^, or equivalently / and E, are varied. Observe also that, from 
the point of view of integrable field theory, the duality M ^ M sends /3^ to l//?'^. As 
M passes through into the 'dual' regime — 1 < M < 0, /3^ increases through 1 and 
the associated perturbation of a c = 1 conformal field theory becomes irrelevant. The 
full range — 1 < M < 00 is mapped onto < /5^ < 00. For M < — 1, the nature of the 
Schrodinger problem ( |2.1D changes fundamentally, as the singularity at the origin ceases 
to be regular. Formally M < — 1 corresponds to purely imaginary values of /?, and it is 
tempting to suppose that there should be a relation with the non-compact sinh-Gordon 
model. However at this stage this is pure speculation. 

9 Conclusions 

In this paper we have further explored the correspondence between certain integrable 
quantum field theories and Schrodinger equations. The two natural classes of spectral 
problems for these equations have been associated with the two types of operators which 
arise in the integrable quantum field theories, namely the T and Q-operators. The 
neatness of this fit convinces us that the matching is no accident, and that it promises 
to lead to a fruitful interaction between these previously distinct fields. A new angle on 
the theory of integrable models seems to be opening up, and this alone should justify 
some further work on the many unresolved issues which remain. 

On a technical level, while we have tried to give most of the details of our calcula- 
tions, our discussion has not been completely rigorous and some aspects will probably 
repay a more careful treatment. In addition, the treatment of this paper should be 
completed with a full discussion of the 'quantum' (non-semiclassical) domain M < 1. 
This is particularly important because the region includes the unitary minimal mod- 
els. One example, the 'Airy case' (M = 1/2, / = 0), has already been treated in j^, 
and its consistency with integrable field theory predictions has been confirmed both 
numerically ||l[ and analytically [^. 

The literature on the spectral theory of the (radial) Schrodinger equation is exten- 
sive, much of it dating back to Regge pole theory, and we can hope to apply further 
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results to the study of integrable quantum field theories. In the other direction, the 
realisation that the some of the functional relations uncovered by Sibuya, Voros and 
others in connection with ordinary differential equations also arise in integrable models, 
and can be solved by means of nonlinear integral equations, promises to be very useful. 
It should be remarked that functional relations can be found for more general potentials 
than the x^^ + /(/+l)x~^ cases used above (see |1J, 2J, 42|). It would be interesting 



to know if these also have a role to play in the theory of integrable models. A relation 
between certain functional relations and differential equations was recently discussed 



in |43], though the approach is rather different from that adopted here. 

The match that we have found connects spectral determinants of the Schrodinger 
equation with the vacuum expectation values of the T and Q-operators. However 
this leaves to one side all of the other expectation values, which emerge when excited 
states in the integrable quantum field theories are discussed. They obey the same sets of 
functional relations as the vacuum expectation values || ,^,27,2^, and a major challenge 
is to find an interpretation for them within the context of Schrodinger equations. One 
way to access excited states in integrable quantum field theories is by a process of 



analytic continuation in a suitable parameter [19|. However if we remain with the 
massless models discussed in this paper, this idea is unlikely to help since the functions 
Tp^ and y are single- valued in the one remaining parameter, namely / (see [14,30]). If the 



correspondence could be extended to cover massive models, then the story would change 
since these models are already known to exhibit a pattern of branch points joining the 



ground state to excited states (examples can be found in |19, 44|). Whether massive 
models can be incorporated into the differential equation framework is another question 
that for the moment remains completely open, and any progress on this issue would be 
extremely interesting. 
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Appendix A WKB details 

In this appendix we sketch how to obtain ( p.21[) . The equation describes the asymp- 
totic behaviour of log-D~(-E', /) = W[y{x, E,l),^p~^{x, E,l)]. Note first that D~ can be 
evaluated as 

D-{E,l) = lim \{2l+l)x^y{x,E,l)\ . (A.l) 

x^O L J 

For convenience we set E = —a, so that the problem under discussion becomes 



^^2+^ +V^ + «b = (A-2) 
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and we are interested in the behaviour as \a\ — > +00. The WKB approximation cannot 
be applied directly to this equation, because of the nature of its singularity at x = jsT 
The remedy found in |3£] is to make the variable change x = e^, y{x) 
Writing A = / + 1/2 the equation becomes 



e^/^0Ui 



+ Riz,a,X) U' = 



with R{z, a, A) = e(2^+2)^ + ae^^ + A^. The WKB approximation for <f> is 



(A.3) 



A 



R{z,a,Xy/^ 



exp 



\/R{u,a, A) du 



+ 



B 



R{z,a,Xy/^ 



exp 



\/R{u, a, A) du 



(A.4) 

In contrast to the WKB treatment of the initial equation, this is approximation is good 
for all real z (and hence for x all the way down to zero) since R~^''^{R~^''^)" tends to 
zero uniformly in z as \a\ — > +00, so long as arga 7^ vr (for a discussion, see |45|]). 

The solution that we hope to approximate is subdominant as x — > +cxd, which 
requires A = 0. The value of B must be chosen so as to match the large-x asymptotic 
{ \2.2\j . A little thought shows that the correctly- normalised approximate solution can be 
written as 



{z,a,X) 



1 



i?(z,a,A)i/4 



exp 



[VRiu,a,X) 



,(M+l)u 



]du 



1 



M+1 



=(A/+l)^ 



(A.5) 



To estimate (A.l) we need to analyse this quantity as z — > —00. An integration by parts 
extracts a term —Xz in this limit, after which A can be dropped from the remaining 
integral, giving a leading a-dependence 



(z,a,A) 



A-1/2 



exp 



A-1/2 



exp 



Xz+ [^e(2A^+2)«+ae2« - e^^^+i)"] du 

POO 

.Xz + a(*^+i)/2^'^ / [Vt^^^ - t^^] dt 
Jo 



(A.6) 



Reverting to the original variables and substituting into (A.l), it is easily checked that 
the result quoted in the main text is recovered. 

If M < 1, the above treatment fails, because the integral in ( |A.5[) diverges. The 
remedy is to introduce further regularising terms, by expanding y^R(u,a, A) to higher 
order. For example, if 1/3 < M < 1 the argument of the exponential in ( |A.5| ) should 
be replaced by 



IVi^R^ - e(*^+y" - ^e(i-*-^)«] du - ^, 



.iM+l)z 



2-2M 



,(1-M)z 



and ( [A.(j ) becomes 

(f){z, a. A) ~ A~ ' exp 



_Az + a(^^+i)/2A'^ 



/•CO 

JO 



(A.7) 



(A. 
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However, for l/3<M<lwe have 

2 / [{f^^ + i)V2 - 1^^ - it-*^] dt = -^r(-i - ^)r(i + ^) , (A.9) 



and so the result (2.22|) is preserved - we have simply continued it around the pole 



at M = 1. Despite the fact that the |i?| — > oo asymptotic of D~ is unchanged, the 
move to M < 1 does have an effect on the limiting behaviour of the auxiliary function 
d{E,l) = uj'^'-^^ D~ {uj"^ E , I) / D~ {uj~'^ E , I) used in §^ The correct choice of branches 
when applying the asymptotic (2.21) at the shifted arguments to E changes, with the 



result that the asymptotic logd{E,l) ~ —^ibo{E)^^ now holds in the revised sector 
~W+i ^ ^'^s{E) < if+T • Outside this sector, the leading contributions cancel and so 
all that can be concluded is that the asymptotic is subleading, though from the nonlinear 
integral equation (|2.32|) we expect that it will be 0(1), and in fact equal to (2/+l)7ri. 
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